Classical antiferromagnet on a hyperkagome lattice 
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Motivated by recent experiments on Na^raOs [Y. Okamoto, M. Nohara, H. Aruga-Katori, and 
H. Takagi, arXiv:0705.2821 (unpublished)], we study the classical antiferromagnet on a frustrated 
three-dimensional lattice obtained by selectively removing one of four sites in each tetrahedron of the 
pyrochlore lattice. This "hyperkagome" lattice consists of corner-sharing triangles. We present the 
results of large- Y mean field theory and Monte Carlo computations on O(N) classical spin models. 
It is found that the classical ground states are highly degenerate. Nonetheless a nematic order 
emerges at low temperatures in the Heisenberg model (N = 3) via "order by disorder" , representing 
the dominance of coplanar spin configurations. Implications for ongoing experiments are discussed. 

PACS numbers: 75.10.Hk, 75.50.Ee, 75.40.Cx 



Antiferromagnets on geometrically frustrated lattices 
often possess macroscopically degenerate classical ground 
states that satisfy peculiar local constraints imposed by 
the underlying lattice structure [l[ . Such highly degener- 
ate systems are extremely sensitive to thermal and quan- 
tum fluctuations, and thereby intriguing classical and 
quantum ground states may emerge via "order by disor- 
der" On the other hand, systems may remain disor- 
dered even at zero temperature 0]. These paramagnetic 
states are called spin liquid phases and their classical 
and quantum varieties have been recent subjects of in- 
tensive theoretical and experimental research activities 
Excitement in such spin systems [3| has also led to 
developments in mesoscopics [5| , optical lattices [6[ , and 
quantum coherence and computing [7]. 

Among several examples of two and three-dimensional 
frustrated magnets, the kagome and pyrochlore lattices 
have obtained particular attention because a relatively 
large number of materials with the magnetic ions sit- 
ting on these lattice structures are available Both 
of these lattices are corner-sharing structures of a basic 
unit; the triangle and tetrahedron respectively. Despite 
this similarity, the classical Heisenberg magnet orders on 
the kagome lattice [3 [p while it remains disordered on 
the pyrochlore lattice [lfjj. The nature of the spin-1/2 
quantum Heisenberg magnets on these lattices has not 
been settled and remains an important open problem 




FIG. 1: (color online). The hyperkagome lattice, 
lines show the underlying pyrochlore lattice. 
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11, 12]. On the other hand, spin-1/2 systems are rare on 



these lattices and other degrees of freedom such as lat- 
tice distortions may play an important additional role. 
As a result, direct experimental tests on spin-1/2 quan- 
tum magnets have been difficult to realize. 

In this context, the recent experiments on Na4lr 3 0s 
[l3| may provide an important clue on these issues, albeit 
in a different three-dimensional frustrated lattice. Here 
Ir 4+ carries spin-1/2 as the five d-electrons form a low 
spin state in the t2 g level. The Ir and Na ions together 
occupy the sites of the pyrochlore lattice such that only 
three of the four sites of each tetrahedron are occupied 



by Ir. The resulting lattice of magnetic Ir is a network of 
corner-sharing triangles as shown in Fig. [TJ where each 
triangle is derived from different faces of the tetrahcdra. 
In analogy to the kagome lattice in two-dimensions, it is 
called the hyperkagome lattice. Even though the Curie- 
Weiss temperature is large, 8\y = — 650K, the suscepti- 
bility and specific heat show no sign of magnetic order- 
ing, nor lattice distortion, down to T ~ |6*;y|/200 [l3j |: 
suggesting that it may be a spin liquid down to low tem- 
peratures. 

In this paper, we study the classical antiferromagnet 
on the hyperkagome lattice. Such investigations not only 
reveal the behavior of the antiferromagnet in the classi- 
cal regime, but also provide an important starting ground 
for the understanding of quantum fluctuation effects. We 
first study the large- N limit of the O(N) vector spin 
model at zero temperature and compute the spin-spin 
correlation function in the large-A^ mean field theory 
14 , 1 X 51 ] . It is found that there exist macroscopically de- 
generate ground states. 

Then we perform large-scale Monte Carlo computa- 
tions on the Ising (N = 1) and the Heisenberg (N = 3) 
models. The Heisenberg model (with exchange coupling 
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FIG. 2: Contour plots of the structure factor in the [hhl] 
plane, (a) large-TV theory at zero temperature, (b) Monte 
Carlo simulations at T/J — 1/100 and L = 8. Axes range 
from — 4-7T to 4tt and both plots are at the same resolution. 



J) remains disordered down to quite low temperatures, 
exhibiting very similar spin-spin correlations to those of 
the large- TV model. These correlations show a character- 
istic dipolar structure in the reciprocal space, which can 
be explained by a mapping to a gauge theory [15l OH. On 
the other hand, the spin correlations in the Ising (TV = 1) 
model turn out to be quite different. 

Most interestingly, a first order transition to a long 
range nematic order is observed in the Heisenberg model 
at a finite temperature. As explained below, this nematic 
order emerges via an "order by disorder" effect and rep- 
resents the dominance of coplanar spin configurations. In 
the disordered phase no evidence of magnetic ordering is 
found while our numerical data cannot definitely confirm 
the presence/absence of magnetic order at temperatures 
below the onset of nematic order. We have also inves- 
tigated the effect of an external magnetic field, h, and 
found that a collinear order is chosen when h = 2J in 
analogy to a similar study on the kagome lattice [1 7I | . 

The lattice and local constraints. — The hyperkagome 
lattice, relevant to Na^^Os, can be represented by the 
simple cubic lattice with a twelve-site basis, as shown in 
Fig. [TJ This lattice is also a three-dimensional network 
of corner-sharing triangles. The model for the classical 
nearest-neighbor antiferromagnet on the hyperkagome 
lattice can be written as 



(ij) A 



constant, 



(1) 



where J > 0, represents the sum over the nearest- 
neighbors, and Si = (S}, Sf) are TV-component spins 
of fixed length TV. Sa = X)ieA ^ s ^ ne vectorial sum 
of the spins in each triangle and J2a represents the sum 
over all triangles. For TV > 2, the classical ground state 
satisfies Sa = for every triangle while the constraint on 
the Ising model is Sa = ±1- Thus, from the outset, one 
may expect that the physics of the larger- TV models would 
be different from the Ising case. This is different from the 
antiferromagnet on the pyrochlore lattice where the Ising 
and Heisenberg models satisfy the same constraint. 



Large-N mean field theory. — Following Refs. 1 14 . Il5l . 
we rewrite the Hamiltonian as H = -j . My Si • 
Sj, where My is the interaction matrix that has 
the information about the nearest-neighbor interac- 
tion. The corresponding partition function is given by 
Z = jVcp VX e - s (<M) with the action S(<p, A) = 

[\ M ij<t>i ■ 4>j + T^yO^i ' 4>i - N )] > where fa is an 
TV-component real vector field and Ai the Lagrange mul- 
tiplier for the constraint (pi ■ (pi = TV. 

Now we take the TV — > oo limit and set a uniform 
A, = Ao- The locations i — (l,fi) of spins can be la- 
beled by those of the cubic unit cell I = 1, ...,n c and 
the lattice sites fi = 1, 12 within the unit cell (n c is 
the total number of the unit cells in the lattice). The 
Fourier transform with respect to the positions of the 
unit cells leads to S = X) q IZ^,i/ I^T^q.M ' ^q,^ witn 
= M£" + S^Xq. Here Ao and the eigenvalues, m^ pi 
of the 12 x 12 interaction matrix M^ v are determined by 

the saddle point equation, 12n c = VV Y^ 2 , ^ , * — . 

It is found that the lowest eigenvalue is four-fold de- 
generate and independent of the wavevector. The next 
lowest eigenvalue has a dispersion and becomes the same 
as the lowest eigenvalue only at g = 0. These fea- 
tures are very similar to those in the kagome and py- 
rochlore lattices. These results imply that the spin 
structure of this system is indeed highly frustrated and 
that magnetic order is suppressed. The static spin- 
spin correlation function can be computed via 14, [l5| 



S-q,i/} 



where Uq tflp is a uni- 



tary transformation that diagonalizes the interaction ma- 
trix, M^ u . At zero temperature, the four degenerate 
eigenvalues dominate the behavior of the spin-spin corre- 
lation function. The resulting zero temperature structure 
factor, 5(q) = ^^(Sq^-S-q^), in the [hhl] plane of the 
reciprocal space is shown in Fig. [2l The presence of high 
intensity along bow-tie structures is apparent and qual- 
itatively similar to that found on the pyrochlore lattice. 
As discussed below, the structure factor in the large-TV 
limit is very similar to that found by Monte Carlo for the 
Heisenberg model (TV = 3) above the nematic ordering 
transition temperature, but quite different from the Ising 
(TV = 1) case. 

Dipolar spin correlations. — We found that the real 
space spin-spin correlation function at long distances is 
well described by the following dipolar form. 



"Jap 



3(e; • r.y)(ej v, 



(2) 



where ry is a vector connecting sites i and j, and a, (3 = 
x,y, z. The "dipolar vectors" are shown in Fig.[3j 

In analogy to the pyrochlore [l5| . we may understand 
the spin correlations in this system by mapping to a pure 
Maxwellian action with a "Gauss law" constraint. We 
first consider a dual lattice of the hyperkagome lattice; 
the sites on the hyperkagome lattice should be placed on 
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FIG. 3: (color online). 
(e' re ) vectors. 



The dipolar (e K ) and dual lattice 



the bonds of the dual lattice. This dual lattice can be ob- 
tained by connecting centers of the tetrahedra of the un- 
derlying pyrochlore lattice, but only along the directions 
that would pass through the sites of the hyperkagome 
lattice. Then there exists a unique bond k (on the dual 
lattice) for a given site i. Let us define e' K as the unit 
vector along the bond k (see Fig. \S§ . We now define N 
number of "magnetic" vector fields b" along the bond 
k via b" = 5*e K , where 5" and e K represent the spin 
and unit vector defined earlier on the hyperkagome lat- 
tice. Notice that the direction of these "magnetic" vector 
fields is not along e' K . Nonetheless the "magnetic" vector 
fields satisfy Div h a = £ K e' K • b« = £ K e' K • e K S£ = 
on the dual lattice. This is a direct consequence of the 
constraint ^2 ieA Sf = and e' K ■ e K = y/2/3 for all k. 

We now define coarse grained "magnetic" vector fields, 
B a , averaged over clusters of spins [la, llo|- There exist 
many "flippable" spin configurations where local rear- 
rangement of the spins in a cluster can be made without 
violating the constraints. The coarse-grained field over 
such "flippable" spin configurations will average out to 
a small value. Then the large entropic weight is related 
to the small values of B a . This feature can be repre- 
sented by an entropic weight of the form exp(— -|- J dr"B 2 ) 
(TEl Hij . This Maxwellian form of the "action" and the 
"Gauss law" constraint will lead to the dipolar form of 
(B? (r)S^(O)) oc 5 aP (ZxiXj - r 2 6 %3 )/r 5 , and hence the 
spin-spin correlation function in Eq. O The discovery of 
such spin correlations supports the entropic argument a 
posteriori. 

Monte Carlo simulations. — Classical Monte Carlo sim- 
ulations for the Heisenberg (N — 3) and Ising (JV = 1) 
models are performed on L x L x L clusters of unit cells. 
We mostly discuss the results of the Heisenberg model 
here and mention those of the Ising model only as neces- 
sary. In the first place, a first-order transition occurs in 
the Heisenberg model. This can be most clearly seen in 
the nematic correlation function defined as |8fl 



(3) 



and n„ 



3^(Si x S 2 + S 2 x S 3 + S 3 x Sr), where Si, 
S2, and S3 are three spins on the triangle a. g(r) = 1 in 
an ideal coplanar state and g(r) = in a non-coplanar 



0.8 
I 0.6 

- 0.4 - 
0.2 
- 
(b) 1 
0.8 
I 0.6 - 
~ 0.4 
0.2 - 
- 



OOO 



oooo^ 



L=9 down 
L=9up 



J 




0.02 0.04 0.06 0.08 

(777) 1/2 



0.1 



FIG. 4: (color online). The nematic correlation function 
for the next-nearest neighbor triangles, (a) Hysteresis is ob- 
served upon lowering (down) or raising (up) the temperature; 
(b) finite size effects scale to a finite first order transition 
temperature (lines guide the eye). 
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FIG. 5: (color online). The energy E(T) - E(T = 0) per spin 
divided by temperature. In the T — > limit, this quantity 
approaches 11/12 and equals the specific heat. 



state. The nematic correlation function for the next- 
nearest neighbors is shown in Fig. [4] and it clearly shows 
a first order transition from a low temperature nematic 
ordered state to a disordered state. Hysteresis associ- 
ated to this transition occurs in the temperature window 
(1 — 5) x 10~ 3 J; coplanar configurations are chosen below 
this window via "order by disorder" . Similar behavior 
is seen for the nearest-neighbor and all higher-neighbor 
correlations. The energy and specific heat data are also 
consistent with the first order transition to nematic or- 
der (see Fig. [5]). Notice that the Monte Carlo data for 
the three largest system sizes L — 6, 8 and 9 are almost 
identical. The zero temperature specific heat approaches 
11/12 per spin which is consistent with the expectation 
that the low temperature phase is dominated by coplanar 
spin structures. Analysis about coplanar states tells us 
that there are 4 quartic and 20 quadratic modes per unit 
cell. Since each quartic (quadratic) mode contributes 
1/4(1/2) to the specific heat Q, the total specific heat 
becomes 11/12 per spin. Interestingly the same zero tem- 
perature specific heat was obtained in the Heisenberg 
model on the kagome lattice Q. The crucial difference 
between two cases, however, is that the nematic order on 
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FIG. 6: (color online). The inverse susceptibility as a func- 
tion of temperature, fit to a Curie- Weiss form for T/J — 
0.9-2.0. 



the hyperkagome lattice is long-ranged at finite temper- 
atures while it becomes long-ranged only in the T — > 
limit on the kagome lattice [8|. 

On the other hand, no magnetic order is seen prior 
to nematic order as we do not find any elastic peaks in 
the spin structure factor. We cannot, however, reliably 
comment whether there is a magnetic ordering or not at 
still lower temperatures. 

The spin correlations in the Heisenberg model in the 
disordered phase are very similar to those in the large-N 
mean field theory (see Fig. [2]) and can be fitted to the 
dipolar form in Eq. O On the other hand, the results 
on the Ising model are markedly different; the spin-spin 
correlation function decays exponentially. We expect this 
is due to the different form of the local constraint in the 
Ising model. 

Magnetization and susceptibility in an external field. — 
The magnetization as a function of an external magnetic 
field, h, is computed by the Monte Carlo simulations on 
the Heisenberg model. At finite temperatures, a weak 
plateau develops at h/J = 2, which leads to a singular 
structure in the susceptibility (not shown). This can be 
explained by the occurrence of a collinear order (up-up- 
down spin structure) by disorder at h/ J = 2 in analogy 
to the kagome lattice case 17| . 

Implications for experiments. — A Curie- Weiss fit to the 
Monte Carlo susceptibility data for T > J, leads to 
Sew — —2.303(5) J (see Fig. [6]). Comparing this with 
the experimental value 9cw — — 650-fT [13j , one obtains 
J k, This suggests that the nematic transition 

may occur around 0.3 — 1.5K if our results are taken 
seriously, and below this temperature coplanar spin con- 
figurations would be preferred. Even though our classical 
computations may not be directly applicable at such low 
temperatures, we suspect that coplanar spin configura- 
tions may still dominate at low temperatures even in the 
quantum regime [l8j]. 

Our results also suggest that the spin correlations at 
T > J/100 ~ 2-3K may be dominated by the physics 
of the classical spin liquid with dipolar spin correlations; 
this will be checked by neutron scattering experiments. 
Notice that there is no sign of magnetic ordering down 



to 2-3 K in the experiment [13J; this may also be consis- 
tent with our Monte Carlo results that show no evidence 
of magnetic ordering above the nematic transition. It 
remains to be seen whether the system develops mag- 
netic order at very low temperatures by "order by dis- 
order" due to classical/quantum fluctuations or prefers 
a magnetically-disordered quantum spin liquid. Finally, 
in the current work, we have not considered the orbital 
degree of freedom which may play an important role in 
the real material. Studies of quantum spin models and 
the role of orbital degrees of freedom, therefore, are im- 
portant subjects of future study. 
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